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Understanding Parton Distributions from Lattice QCD: 

Present Limitations and Future Promise * 
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This talk will explain how ground state matrix elements specifying moments of quark 
density and spin distributions in the nucleon have been calculated in full QCD, show how 
physical extrapolation to the chiral limit including the physics of the pion cloud resolves 
previous apparent conflicts with experiment, and describe the computational resources 
required for a definitive comparison with experiment. 

1. INTRODUCTION 

The focus of this workshop, the QCD structure of the nucleon, has two complementary 
aspects. One is measuring the quark and gluon structure experimentally. Decades of 
impressive high energy scattering experiments utilizing our knowledge of perturbative 
QCD and factorization have produced detailed knowledge of the light cone distributions of 
quarks and gluons in the nucleon, with the promise of even richer phenomenology to come. 
The second, and presently less well developed aspect, is calculating and understanding 
this quark and gluon structure theoretically, and the only method known at present to 
solve nonperturbative QCD is lattice field theory. Both aspects are essential, and I will 
show that the tools for quantitative calculation of nucleon structure are at hand, explain 
the scale of computational resources required for a definitive calculation, and argue that 
our field needs to provide these computational resources as well as experimental resources 
if we are to truly understand hadron structure. 

Using the operator product expansion and lattice field theory, it is possible to calcu¬ 
late moments of quark distributions, and I will discuss here the first calculations in full 
qcd mg. A major puzzle in the field has been the fact that quenched calculations of 
these moments, which ignore quark-antiquark excitations of the Dirac sea, disagree with 
experiment at the 20-50% level. I will show that contrary to earlier conjectures, at the 
quark masses accessible in practical calculations, including quark loops does not alter 
the results significantly. Rather, I will argue that the physical origin of the discrepancy 
with experiment has been incorrect extrapolation to the physical quark mass, and will 
show how extrapolation incorporating the leading non-analytic behavior required by chiral 
symmetry produces consistent results for the moments of quark distributions. In addi¬ 
tion, we have also compared full QCD results with configurations that have been cooled 
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Figure 1. Connected (upper row) and disconnected (lower row) diagrams contributing to 
hadron matrix elements. The left column shows typical contributions of quarks and the 
right column shows contributions of antiquarks 

to remove all the gluon contributions except for those of instantons and shown that the 
qualitative behavior of the moments is reproduced by the instanton content of the gluon 
configurations. 

2. MOMENTS OF QUARK DISTRIBUTIONS IN THE PROTON 

By the operator product expansion, moments of the following linear combinations of 
quark and antiquark distributions in the proton 

(x n ) q = j dxx n (q(x) + (~l) n+1 q(x)) (1) 

J 0 

{x n )A q = [ dxx n (Aq(x) + (-l) n Aq(x)) 

Jo 

{x n )s q = I dxx n (Sq(x) + (~l) n+1 5q(x)) 

Jo 

are related to the following matrix elements of twist-2 operators 

(PS^^HD^ ■■■iD^^\PS) = 2(x n ~ 1 ) q P {lil (2) 

(PS|^ 7 {mi 7 5 UU 2 • • • iD^^PS) = 2(x n - 1 ) Aq MS {til P^ ■ ■ ■ P^ } 
{PS\^a [a{ ^^ b iD^ ■ ■ ■ xD^ } ^\PS) = 2(x n - 1 )s q MS [a P {fll] P^ ■ ■ ■ P^ } . 

Here, q = q^ + q±, A q — q^ — q^, Sq — qr + q±, x denotes the momentum fraction carried by 
the quark, and { } and [ ] denote symmetrization and antisymmetrization respectively. We 
note that odd moments ( x n ) q are obtained from deep inelastic electron or muon scattering 
structure functions F 1 or P 2 , even moments of (x n ) Aq are determined from g 1: and these 
moments are proportional to the quantities v n+ \ and a n defined in ref. |3j . In addition, g 2 
also determines the quantity d n 

(PSl^^iD^ • --iD^MPS) = -d n S [(T P { ^ ] • • -P^ } 

n 


(3) 
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Figure 2. Plateaus in measurements of operators in a zero momentum ground state as a 
function of the euclidean time separation from the source. 

although the lattice calculation of d n with Wilson fermions is complicated by mixing with 
the lower dimension operator “7^757^ ' •' iD^ n \ 

Even moments ( x n ) q are obtained from deep inelastic neutrino scattering, and in ad¬ 
dition, a variety of other processes have contributed to what is now a detailed empirical 
knowledge of the quark and antiquark distributions in the nucleon. We will subsequently 
compare our results with moments calculated from the CTEQ, MRS, and GRV global fits 
to the world supply of data. 


3. CALCULATION OF MATRIX ELEMENTS 

Proton matrix elements of the operators in Eq.^j are calculated by evaluating the con¬ 
nected and disconnected diagrams shown in Fig. |l|. Note that both the connected and 
disconnected diagrams each receive contributions from quarks and antiquarks. Depending 
on the moment, by Eq.|l], the sum of the diagrams yields either the sum or difference of the 
moments of the quark and antiquark distributions. Because it is technically much more 
difficult to evaluate the disconnected diagrams, our present calculations only include con¬ 
nected diagrams. Fortunately, the disconnected diagrams are flavor independent, so they 
cancel out of the flavor non-singlet difference between up and down quark distributions. 
Hence, in Table we compare lattice calculation of the difference between connected 
diagrams for up and down quarks with the corresponding difference of moments of exper¬ 
imental data for the sum or difference of quark and antiquark distributions. 

On the lattice, connected diagrams are evaluated by calculating a three point function 
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Table 1 

Lattice Operators 



H(4) mix 

P 

lattice operator 

(x)^ 

63 

no 

7^0 

?7{i Dx\q 

( X )q ] 

3 1 

no 

0 

| E-=i QliDiq 

(x 2 )q 

sr 

yes 

7^0 

q^liDiDijq - 5 Ef = 2 h'p/L/Lp/ 

(x 3 ) q 

2 1 

no* 

7^0 

^{iDiDAD^q + q r )'{2D2D3D3}q ~ (3 4) 

< 1 )a. 

4 i 

no 

0 

<Z7 5 73<? 

( X )aI 

63 

no 

7^0 

q^ 5 ^{iD 3 }q 

( x )aI 

63 

no 

0 

<77 5 7{3.D4}g 

( x2 )a q 


no 

7^0 

g7 5 7{iD 3 D4}g 

(!)Sq 

6 1 

no 

0 

g7 5 cr 34 g 

( X )sq 

sr 

no 

7^0 

q^^Diyq 

d\ 

6 1 

no** 

0 

g7 5 7[sD4]g^ 

G ?2 

sr 

no** 

7^0 

g7 5 7[iD { 3iP>4}g 


in which a source creates a state with the quantum numbers of the proton, the operator 
acts on this state, and a sink finally annihilates the state. Because evolution in imaginary 
time Liters out the ground state, when the operator is sufficiently far from both the source 
and sink, it acts in the ground state and produces the desired ground state matrix element. 
As the time at which it acts approaches either the source or sink, it sees excited state 
contaminants, yielding a central plateau corresponding to the physical matrix element 
and exponential contaminants at the edges. Obviously, it is beneficial to optimize the 
overlap of the source with the ground state to maximize the plateau region and minimize 
the effect of the excited state contaminants at the edges. 

In this work, connected diagrams were calculated using sequential propagators gener¬ 
ated by the upper two components of the nucleon source J a = e abc . The 

overlap with the physical proton ground state was optimized using Wuppertal smearing 
[0] to maximize the overlap P(0) = |(J|0)| 2 . Varying the smearing reduced P by over 4 
orders of magnitude, yielding an overlap with the physical ground state of approximately 
50%. Dirichlet boundary conditions were used for quarks in the t-direction. 

The resulting plateaus for four operators that could be measured in a proton with 
zero three-momentum are shown in Fig. |0. Here one observes both a statistically well 
determined central plateau region and the effectiveness with which the excited state con¬ 
taminants have been reduced by the optimized source. 

4. OPERATORS AND PERTURBATIVE RENORMALIZATION 


The continuum operators defined above are approximated on a discrete cartesian lattice 
using representations of the hypercubic group that have been chosen to eliminate operator 
mixing as much as possible and to minimize the number of non-zero components of the 
nucleon momentum. The operators we have used are shown in Table 1, where we have 
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Table 2 

Perturbative renormalization 


jdLATT jdMS 7 7 

jD_£>_ ^/3= 6.0 -^/3=5.6 

>1U/I8fi _ 40 n 080 0 088 


<■ X >i aj 

8 

3 

-3.16486 

40 

9 

0.989 

0.988 

w 

8 

3 

-1.88259 

_40 

9 

0.978 

0.977 

(X 2 ) q 

25 

-19.57184 

— §1 

9 

1.102 

1.110 

EE 

157 

30 

-35.35192 

2216 

225 

1.215 

1.231 

(1)a, 

0 

15.79628 

0 

0.867 

0.857 

(x)aI 

8 

3 

-4.09933 

40 

9 

0.997 

0.997 

eS 

8 

3 

-4.09933 

40 

Sl 

0.997 

0.997 

E)a q 

25 

6 

-19.56159 

_ 67 

9 

1.102 

1.110 

(X)iq 

l 

16.01808 

-l 

0.856 

0.846 

( X )sq 

3 

-4.47754 

-5 

0.996 

0.995 

d\ 

0 

0.36500 

0 

0.997 

0.997 

d 2 

7 

6 

-15.67745 

35 

ia 

1.116 

1.124 


indicated whether the spatial momentum components are non-zero and whether mixing 
occurs. Note, no* indicates a case in which mixing could exist in general but vanishes 
perturbatively for Wilson or overlap fermions and no** indicates perturbative mixing 
with lower dimension operators for Wilson fermions but no mixing for overlap fermions. 
Because the statistical errors are much larger for sources projected to non-zero momentum, 
the moments corresponding to operators requiring non-zero momentum are presently not 
well determined. 

To convert from lattice regularization at the scale of the inverse lattice spacing 1/a 
to the continuum MS scheme at momentum scale Q, we use the one-loop perturbation 
theory result 


or s m =e ( e - + 
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The anomalous dimensions and the finite constants Bij we have calculated and used 
in this work are tabulated in Table 2 || . 


5. RESULTS 

The moments listed in Table 1 were calculated |I],|2j on 16 3 x 32 lattices for Wilson 
fermions in full QCD at [3 — 5.6 using 200 SESAM configurations at each of 4 k's and 
at f3 = 5.5 using 100 SCRI configurations at 3 k's. They were also calculated with two 
sets of 100 full QCD configurations cooled with 50 cooling steps and in quenched QCD 
at (3 — 6.0 using 200 configurations at each of 3 k's. Typical linear chiral extrapolations 
for operators calculated with nucleon momentum equal to zero are shown in Figure § for 
full and quenched calculations of (x) q and (x)a q , showing agreement within statistical 
errors. To avoid finite volume errors at the lightest quark mass, the SESAM || results 
were extrapolated using the three heaviest quark masses. Table 3 shows a major result of 
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Figure 3. Comparison of linear chiral extrapolations of full and quenched calculations of 
(x)^ and (1 )aij showing agreement within statistical errors 


our work, that there is complete agreement within statistics between full and quenched 
results. Statistics with the SCRI configurations |7|] are currently being increased to study 
the dependence on the coupling constant. 

Typical chiral extrapolations for cooled configurations are compared with the corre¬ 
sponding uncooled full QCD calculations in Figure f|. This qualitative agreement be¬ 
tween cooled and uncooled results occurs at light quark mass for all the twist-2 matrix 
elements we calculated and demonstrates the degree to which the instanton content of 
the configurations and their associated zero modes dominate light hadron structure |5|. 

A major crisis and longstanding puzzle in this field has been the fact that when quenched 
lattice results are linearly extrapolated in m qi results disagree with experiment at the 20- 
50% level. Our results clearly show that inclusion of quark-antiquark excitation from the 
Dirac Sea does not resolve this discrepancy. Salient examples from Table |3| are ( x) u -d ~ 
0.25 — 0.27 compared with the experimental result 0.15 and qa = (l)Au-Ad ~ 1-0 — 1.1 
compared with the experimental result 1.26. 




Figure 4. Comparison of linear chiral extrapolations of full and cooled calculations of 
(x)^ and (I)a,} showing the extent to which instantons reproduce the full result. 
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Table 3 

Comparison of linear extrapolations of our full QCD and quenched results with other 


lattice calculations and phenomenology at 4 GeV in MS 


Conn. 

M. E. 

QCDSF 

QCDSF 
(a = 0) 

SESAM 

Quenched 

Full QCD 
(3 pts) 

Phen. 

(q±q) 

{x) u 

0.452(26) 



0.454(29) 

0.459(29) 


{x)d 

0.189(12) 



0.203(14) 

0.190(17) 


(%)u—d 

0.263(17) 



0.251(18) 

0.269(23) 

0.154(3) 

( X 2 ) u 

0.104(20) 



0.119(61) 

0.176(63) 


(x 2 ) d 

0.037(10) 



0.029(32) 

0.031(30) 


{X 2 ) u -d 

0.067(22) 



0.090(68) 

0.145(69) 

0.055(1) 

( X 3 )u 

0.022(11) 



0.037(36) 

0.069(39) 


(x 3 ) d 

-0.001(7) 



0.009(18) 

-0.010(15) 


( X 3 ) u -d 

0.023(13) 



0.028(49) 

0.078(41) 

0.023(1) 

(1 )a« 

0.830(70) 

0.889(29) 

0.816(20) 

0.888(80) 

0.860(69) 


(1 )a d 

-0.244(22) 

-0.236(27) 

-0.237(9) 

-0.241(58) 

-0.171(43) 


(l) Au—Ad 

1.074(90) 

1.14(3) 

1.053(27) 

1.129(98) 

1.031(81) 

1.248(2) 

{x)au 

0.198(8) 



0.215(25) 

0.242(22) 


{x)a d 

-0.048(3) 



-0.054(16) 

-0.029(13) 


{x^j Au—Ad 

0.246(9) 



0.269(29) 

0.271(25) 

0.196(9) 

(x 2 )a u 

0.04(2) 



0.027(60) 

0.116(42) 


{x 2 )a d 

-0.012(6) 



-0.003(25) 

0.001(25) 


{x ) Au—Ad 

0.05(2) 



0.030(65) 

0.115(49) 

0.061(6) 

Su c 

0.93(3) 

0.980(30) 


1.01(8) 

0.963(59) 


5d c 

-0.20(2) 

-0.234(17) 


-0.20(5) 

-0.202(36) 


d u 2 

-0.206(18) 



-0.233(86) 

-0.228(81) 


dt 

-0.035(6) 



0.040(31) 

0.077(31) 



A breakthrough in resolving this crisis has been the understanding that the physical 
origin of these discrepancies is the inadequate treatment of the pion cloud in the nucleon 
that has been necessary because of limited computational resources. By necessity, present 
calculations are restricted to quark masses that are so heavy that the pion mass is above 
600 MeV and the pion cloud surrounding the nucleon is strongly suppressed. Physical 
quantities like the nucleon magnetic moment and axial charge clearly depend strongly 
on the pion current, and should therefore be very sensitive to the absence of the full 
pion cloud. Furthermore, because of the rapid, nonlinear variation from the chiral logs 
arising from Goldstone boson loops, it is clear that a linear extrapolation is completely 
inadequate to describe the correct chiral physics. 

we have shown that chiral extrapolation incorporating the lead- 


in a recent work [10 


ing non-analytic behavior from chiral perturbation theory can systematically resolve the 
discrepancies in the moments (. x n ) u -d using the formula: 


+ b n ml (4) 

The coefficient of the leading non-analytic behavior m^ln(m^) is determined from chiral 
perturbation theory. The parameter fi specifies the scale above which pion loops no longer 
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Figure 5. Chiral extrapolation of (x)^ , ( x 2 ) q , and ( x 3 ) q using Eq. [| 


produce rapid variation. It corresponds to the upper limit of the momentum integration 
generated physically by the inverse size of the quark core of the nucleon that serves as 
the source for the pion field. As shown in Fig. [|, the value // ~ 550 MeV from ref [|l(] 


which is consistent with the value required to extrapolate the nucleon magnetic moment 
and with chiral nucleon models, extrapolates the world’s supply of lattice data to the 
experimental values of (x) u _ d , {x 2 ) u -d , and (x 3 ) u _d- From this figure, it is also clear that 
5% measurements down to m 2 = 0.05 GeV 2 , requiring a spatial volume of order (4.3 fm) 3 , 
would provide data for a definitive lattice calculation. From the known scaling of present 
hybrid Monte Carlo algorithms, this calculation will require 8 Teraflops-years and thus 
can be carried out on the next generation of 10-Teraflops computers. 


FUTURE PROMISE 

With the present tools of lattice field theory and the physical understanding of the 
essential role of the pion cloud, we now have the framework for definitive calculation of 
hadron structure. Improved actions will allow a reliable extrapolation to the continuum 
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Figure 6. Graph showing how the propagation of a pion in the QCD vacuum (diamonds) 
is well approximated by a propagator containing only the lowest 128 eigenmodes 


limit. Chiral fermions now implement an exact chiral symmetry on the lattice. Partially 
quenched chiral perturbation theory provides a systematic framework to measure the rel¬ 
evant parameters of effective field theory for unambiguous chiral extrapolation using opti¬ 
mal combinations of valence and sea quark masses, including corrections for finite volume 
and fixed topology. As explained above, computers providing sustained performance of 
tens of Teraflops will enable calculation in the chiral regime of sufficiently low quark mass 
and sufficiently large volume for reliable chiral perturbation theory extrapolation. Renor¬ 
malization coefficients and mixing parameters can now be calculated nonperturbatively. 
And finally, eigenmode expansion techniques jR1-|TT| using the dominance of propagators 
by low eigenmodes as illustrated in Fig. |6] enable the calculation of previously intractable 
disconnected diagrams. 

In conclusion, understanding the quark and gluon structure of QCD requires not only 
frontier facilities for experimental measurements, but also frontier facilities for fundamen¬ 
tal lattice calculations. The time has come to think of dedicated computers in the same 
way as one thinks of experimental apparatus for the field. And, as has been so effective 
in experimental physics, the time has also come for international collaboration in large 
scale computation. 
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